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We determine the current -voltage characteristic of type II superconductors in the presence of 
strong pinning centers. Focusing on a small density of defects, we derive a generic form for the char- 
acteristic with a linear flux- flow branch shifted by the critical current (excess-current characteristic). 
The details near onset, a hysteretic jump (for k ^ 1) or a smooth velocity turn-on («: — >■ 1), depend 
on the Labusch parameter k characterising the pinning centers. Pushing the single-pin analysis into 
the weak pinning domain, we reproduce the collective pinning results for the critical current. 
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The defining property of a (type II) superconductor is 
its ability to carry electric current without dissipation. 
This superflow is destroyed when the magnetic induction 
B enters the material in the form of quantized flux lines 
or vortices [1]: driven by the current density j via the 
Lorentz force Fl — jB/c^ the finite velocity v of vortices 
generates a dissipating electric field E — vB/c parallel 
to j [2]. It is the material defects immobilizing vortices 
which reestablish the superflow of current, eventually 
rendering the superconductor amenable to technological 
applications. An elementary distinction is made in the 
design and action of pinning defects: strong pins act in- 
dividually and generate large (plastic) deformations and 
metastable vortex states, while weak defects are unable 
to pin vortices alone and thus act collectively. In this let- 
ter, we determine the generic force-velocity (or current- 
voltage) characteristic of vortices driven by a current j 
and subject to a small density of strong pins. 

Vortex pinning has originally been studied by Labusch 
[3] for strong pins (see also Ref. 4) and has later been ex- 
tended to weak collective pinning by Larkin and Ovchin- 
nikov [5]. While the latter has been profoundly studied 
[6, 7], the further development of strong pinning theory 
has been less dynamic, although some progress has been 
made over time [8-12]. Recently, the two regimes have 
been analyzed within a pinning diagram [13] delineating 
the origin of static critical forces Fc as a function of de- 
fect density Up and strength fp. Here, we go beyond the 
calculation of the static critical force and determine 
the full force i^L (= jB/c) versus velocity v {= cE/B) 
(or j~E) characteristic of a so-called 'hard' type II su- 
perconductor. We focus on the single-pin-single-vortex 
strong pinning regime, implying that defects are dilute 
and moderately strong, pinning only one vortex line at a 
time; furthermore, we concentrate on isotropic material 
and ignore effects of thermal fluctuations. 

The calculation of critical forces for weak pins involves 
dimensional [5, 6] or perturbative [14, 15] estimates and 
is rather on a qualitative level. Calculations of the force- 
velocity characteristic focus either on the perturbative 
regime at high velocities [14, 15] or on the universal 
regime near depinning [16]. The situation is different for 



strong pinning: here, the critical force and the full dy- 
namical response can be determined quantitatively, once 
the shape of the pinning potential is known. The force- 
velocity characteristic we find agrees well with numer- 
ous (even textbook [17, 18]) experimental results [19- 
21]: a nearly linear flux-fiow curve shifted by the critical 
force Fc (excess-current characteristic) , with a hysteretic 
jump in velocity at onset for strong pinning changing to 
a smooth rise on approaching the weak pinning domain. 
Quite remarkably, continuing our single-pin analysis into 
the weak pinning domain, we can find the usual weak 
collective pinning results for the critical force. Below, 
we derive the formalism leading us to the force-velocity 
characteristic, present the results for the average pin- 
ning force {Fp)(v) for a Lorentzian-shaped pin, derive the 
generic characteristic for the strong pinning case in the 
dilute-pin limit, and finish with a rederivation of the weak 
collective pinning results for the critical current from a 
study of single-defect pinning. 

The velocity-force characteristic derives from the dy- 
namical equation for vortex motion 

7jv^F,{j)-{Fp){v) (1) 

with the Bardeen-Stephen [2] viscosity rj ^ BHc2/ PnC^ 
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FIG. 1. Average pinning force (fp) for a Lorentzian shaped 
pinning potential of various strengths. For strong pinning 
K 3> 1 the critical force /c is large and the pinning force decays 
monotonously. On approaching the Labusch point k — > 1 
the critical force fc vanishes and the pinning force is non- 
monotonic, first increasing oc ^/v and then decaying cx 1/ y'w. 
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{Hc2 is the upper critical field and p„ denotes the nor- 
mal state resistivity) and the velocity dependent aver- 
age pinning force density {Fp){v) (we choose {Fp){v) to 
be positive). The pinning force density is determined 
by the positions (R^ -|- Ufj,{z,t),z) of the flux lines (we 
choose — (R^, z) to be the static lattice positions and 
u^(z,i) is the vortex displacement field; the vortex den- 
sity is Oq^ — B/^o with $0 = hc/2e the flux quantum) 
and the individual pinning forces fp(R — Ri)(5(z — z^) of 
defects located at positions = (Rj,Zi), 



1 ^ r dz 
Fp> - ^^Ey :^Fp(^m,UaO, with 



(2) 



Fp(r^,u^) = ^ ^fp[Rp + u^{z,t) - Ri]S{z - Zi). 

For a point-like defect — ep5(r), the convolution with the 
vortex shape l-|*(r)p w 2^^ / {R^ + 2(^) (* denotes the 
complex order parameter and ^ is the coherence length) 
provides the pinning potential 

ep(R, z) = -ep-^^-p^(5(z) = ep(R) (5(z), (3) 

which here is of Lorentzian shape but may have another 
form in general. The pinning force is given by the gradi- 
ent fp(R) = -VR,ep(R). 

The calculation of the pinning force density requires 
knowledge of the displacement field u^(z,t). The latter 
is determined by the solution of the dynamical equation 
which we write in integral form 

u^{z,t) =vt + al^ j dz'dt' G{'R^-'Rf,,z- z,t-t') 

xFp[r;,u^(z',i')]- (4) 

The first term accounts for the Lorentz force in Eq. (1) 
generating the flux- flow velocity v = Fi^/tj in the ab- 
sence of pinning. The dynamical elastic Green's function 
G{r,t) is given by the Fourier transform of the matrix 



G'a/3(k,a;) 



irjoj 



(5) 



CgqK^ + Ciiki — irjuj 

with the elastic moduli cn (compression), C44 (tilt), and 
C66 (shear) [6]. The task simplifies considerably in the 
dilute-pin limit (to order Up) and for moderately strong 
pinning defects trapping no more than one vortex; in this 
situation we can drop the sums over i and /i in Eqs. (2) 
and (4) . We choose the pin position at the origin and let 
the vortex move on the x-axis; the problem then reduces 
to the calculation of the displacement field Ux{z — 0,t) 
at z = 0. With the asymptotic position x = vt [dX z = 
±00), we have to solve the self-consistent equation 



u{x) = X + 



^ dx' 



G[Q,{x-x')/v]fp[u{x% (6) 



where u{x) = Ux{z = Q,t) and with G = Gxx and fp 
the force along x. Inserting the solution back into Eq. 
(2), we obtain the pinning force density _Fp(0, z,u) — 
— fp[u(vt)]6{z). The average pinning force density 
Upa^J dz{Fp) due to a finite density Up of defects in- 
volves the average (•) over pin locations and time; the 
latter transform to an average along x and the impact 
parameter b of the vortex on the defect, 



{Fp){v) = npifp) = ~np(^ 



°° dx \ 

-fp[u{x)]) (7) 

an / b 



where a|| is the distance between vortices along the x 
direction. Restricting ourselves to the case of strongest 
pinning with 6 = and treating all trajectories within the 
range u ~ ^ of the pin equally, the average over impact 
parameters {■)b contributes a factor a/a± with a_L the 
transverse distance to the next vortex, hence a\\a± — Oq. 
Inserting the result for {Fp){v) back into the dynamical 
equation (1) and solving for the velocity v for a given 
current j provides us with the desired result, the force- 
velocity characteristic of the superconductor. 

In the static situation, the self-consistent integral equa- 
tion (6) simplifies to the algebraic equation 



Us{x) = X + fp[Us{x)]/C, 



(8) 



with ~ G{r = 0,uj = 0) the local static elastic 

Green's function, C ^ So/qq with Eq = {^q/AttX)'^ the 
energy scale for vortices. Strong pinning is character- 
ized by the appearance of bistable solutions in Eq. (8), 
implying that the derivative 



dx 



fp[Us{x)] = -C 



(9) 



of the effective force has to diverge — this provides us with 
the Labusch criterion [3] k = max^ {fp[us{x)]/C} = 1 
separating weak (k < 1) and strong (k > 1) pinning. 
Note that the effective force gradient inside a very strong 
pin is universally given by the effective elastic constant 
C, not by fp. The different solutions of Eq. (8) at k > 1 
are associated with the unpinned {us outside the pin) 
and pinned {ug inside the pin) states of the vortex; their 
asymmetric statistical occupation at finite drive produces 
a finite critical force density F^ = max{Fp) {v — 0) where 
the maximum is taken over the pinned and unpinned 
branches. For a weak pin (k < 1), Eq. (8) has a unique 
solution and the critical force density Fc vanishes. 

In the dynamical situation with v > we have to solve 
the self-consistent integral equation (6) and thus need 
to know the time dependence of the Green's function 
G{0,t). At short times t < tth = V^'o/^^cqq the Green's 
function is dominated by the response of an individual 
vortex line (the one-dimensional (ID) regime), G^°(0, t > 
0) ^ {tth/tY''^ /Cttw, at long times t > ithAVoo the full 
3D vortex system provides the response and G'^°(0,t > 
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0) ~ (tth/i)^^^ao/C'ithA (the intermediate dispersive or 
4D regime with ^ {tth/t)"^ /Ctth is less relevant in 
our analysis below). 

At high velocities v the time integral in Eq. (6) extends 
over short times and the velocity-dependent part of the 
pinning force scales as tG^° oc ^/i, while, at small ve- 
locities, long times are relevant and tG^° cx \/l/t. The 
time t to velocity v transformation t ^ Ucs/v involves 
the effective pin size (Jeff ~ '^^/(l + '^/i^'^th) which de- 
pends on the pinning strength k, and on the velocity v 
itself [22]; for k — > 1 and at high velocities v > K^Vth 
the effective pin size saturates at the true geometric 
pin size f (here, Wth ~ f/^th is the basic velocity scale). 
The corrections to the critical force Fc at small velocities 
V < (ao/A^)Ktith then are expected to scale as \/v/KVth, 
while the high-velocity v > K^Vth corrections to the dissi- 
pative force rjv (flux-flow) decay as \fv^^Jv. This is con- 
firmed by the numerical solution of the problem follow- 
ing the steps indicated above and where the results are 
shown in Fig. f (we assume non-dispersive moduli corre- 
sponding to a field B ~ $o/A^). The forward integration 
of Eq. (6) has been done for a Lorentzian-shaped poten- 
tial of the form (3) and different pinning strengths as 
expressed by the Labusch parameter k ^ (ep/f£o)(ao/0; 
with Bp ~ ^ £oC {He the thermodynamic critical 

field) the Labusch parameter can naturally access large 
numbers k ^ oq/^ 3> f . The scaled average pinning force 
{fp)a^/ep£_ is plotted against the scaled velocity v/nvt^ 
and exhibits a monotonic decrease at large k, and a non- 
monotonic behavior enforced by the vanishing of the crit- 
ical force fc — {fp){0) as K — > 1. While our rough esti- 
mate above correctly predicts the shape oc \/v/ KVth of 
the finite-velocity corrections, its sign depends on k in a 
nontrivial way [22]. Note that we plot the single-pin re- 




FIG. 2. Force-velocity characteristic for a Lorentzian-shaped 
pinning potential. For small defect densities ripao^^K — 0.05 
we find an excess-current characteristic, a shifted (by Fc oc 
Up) linear curve with a slope reflecting fiux-flow behavior and 
approaching the true (unshifted) fiux-flow behavior only at 
high velocities nvth ^ Vc- The insets sketch the behavior near 
critical, with a hysteretic jump of order v^i oc appearing 
at strong pinning k ^ 1 and a smooth onset v oc Vth{n — 
1)* (Fl/Fc — 1)^ on approaching the Labusch point k — )> 1. 



suit (/p) rather than the corresponding force density (Fp) 
as the density Up is an important independent parameter. 

In our discussion of the force-velocity characteristic we 
first concentrate on the overall shape away from the onset 
of vortex motion. The generic characteristic 

^ ^ ^ _^ {fp)iv/liVth) ^^Q^ 
Fc Vc fc 

involves two velocity scales, the velocity KVth govern- 
ing the pinning force (fp) (as confirmed by a detailed 
analysis of Eq. (7) [22]) and the scale Vc — Fc/rj ap- 
pearing from the competition between the dissipative 
{ijv) and the critical {Fc — Upfc) force densities. In 
the limit of small pin densities rip, the linear term in 
Eq. (10) changes on the small velocity scale Vc oc Up, 
while the pinning force {fp){v) deviates from its static 
value fc only on the larger scale KWth which does not 
depend on rip. Indeed, making use of the expression 
Pc ^ {S.^ /o.D'iT'pfpiK — 1)^/k for the critical force density 
[13], we find the ratio Vc/nvth ^ npao^'^iK—l)'^/K ^ I in 
the small pin density limit and at fixed k (note that the 
limit K — ^ I at fixed rip would take us out of the single-pin 
regime). With (/p)(w) « fc for velocities v ^ Vc <^ nvth 
we find a characteristic that takes the generic form of a 
shifted (by Fc) linear (flux-flow) curve, v « (Fl — Fc)/ri, 
see Fig. 2; the free dissipative flow v — Fi^/rj is ap- 
proached only at very high velocities v ^ KWth S> Wc- 
The simple excess-current characteristic is a consequence 
of the separation of velocity scales Vc and KWth; the latter 
merge at strong pinning with increasing density Up when 
strong 3D pinning goes over into ID strong pinning at 
UpttQ^'^K ~ I [13]. Using qualitative arguments, a similar 
excess current characteristic has been found in Ref. 8. 

The above simple overall structure of the force-velocity 
characteristic is modified at very small velocities and 
in close vicinity to the critical force density Fc] in 
this regime we can rewrite Eq. (10) in the simple form 
F^/Fc = v/vc + 1 ± {v/v^f/^, where the '-)-' ('-') 
sign applies to the limits n ^ 1 (k ^ 1). The 
small-velocity pinning scales derive from the 3D ex- 
pression of the pinning force density [22] (Fp) — Fc ~ 
(^2/ao A)np/pKA/ v/ KVth (at large k), ^ {\^/al)KVth 
for K > 1 and u+ ~ (A^/aQ)(K — l^Vth for k — >■ 1. For 
strong pinning, the negative (non-linear) correction in the 
average pinning force density generates a bistability (and 
hence hysteretic jumps) on the scale v-n\ ~ v1/v~ oc n^. 
On the other hand, approaching the Labusch point, 
the correction changes sign and the velocity increases 
quadratically v ~ v^{Fi^/Fc — 1)^ until crossing over into 
the linear regime at ~ ^l/Vp ^ Vc ^ Vp . These 
features are visible in the insets of Fig. 2 showing an 
expanded view of the characteristic near onset. 

Next, we push our single-pin (SP) analysis into the 
weak pinning domain k < 1 and establish its relation 
to weak collective pinning (WCP) theory. In the dy- 
namical formulation of WCP, we determine the pinning 
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force {Fp)^'^^{v) perturbatively (to lowest order in k 
and Up) at high velocities and follow the velocity cor- 
rection 6v = {Fp)^'^^ {v) I r] down to small v. As the 
correction 5v becomes of order v, higher order terms be- 
come relevant [15] and we stop the analysis, interpreting 
the breakdown of perturbation theory as the signature 
of a finite critical force i^™'-"'. The latter then derives 
from the critical velocity Vc defined through the criterion 
(Fpr^^(z;,)^r,i;e = i^™°". 

Within the SP analysis valid at small densities rip, we 
usually calculate the pinning force density {FpY^[v) ex- 
actly, cf. Fig. 1; in the case of weak pinning k <^ 1, we 
can use perturbation theory as well and we find the result 



{Fpriv) 



dtG{0,t)K''"°'{vt,0), (11) 



where we have expanded the average pinning force den- 
sity Eq. (7) for a displacement u{x) = a; -I- Su{x) close 
to flux-flow and used the lowest order (in k) approx- 
imation of Eq. (6) for Su{x). In Eq. (11), K{u) — 
{nplal) J d^i?ep(R — u) ep(R) replaces the usual pinning 
energy correlator showing up in WCP theory [G] (the su- 
perscripts denote derivatives with respect to Ux and Ua). 
Hence, the corrections Sv from both the WCP- and the 
SP analysis agree with one another to lowest order in 
K and in the pin density rip. The difference in the two 
approaches arises when we take the velocity v to zero: 
While we stop a,t 5v ^ v and arrive at a finite F^'^^ in 
WCP, we take v all the way to zero within the SP analy- 
sis and obtain a vanishing critical force = 0. On the 
other hand, using the SP result Eq. (11) and adopting 
the WCP cutoff, we find a finite critical current as well: 
with the estimate (Fp)^^(r;) - np{^/X){f^/so)iv/vth)^/^ 
valid at low velocities and the conditions {Fp)^^{vc) ~ 
rjVc ^ jcB / c, we obtain the critical current 



(12) 



in agreement with the results obtained from weak collec- 
tive pinning theory [13]. This result is quite remarkable: 
first, the critical current (12) is proportional to rip, the 
square of the pin density rip, i.e., its origin is in the corre- 
lations between pins. Second, the result is still consistent 
with the standard SP result {FpY^{v = 0) = 0, as the 
latter is an order rip result and corrections cx rip are be- 
yond the standard SP approach. Going back to strong 
pinning k > 1, we already obtain a finite critical force 
{FpY'^{v = 0) cx rip, linear in pin density. Pin-pin cor- 
relations then are expected to provide corrections o(np) 
which vanish faster than linear and we can approach the 
critical force parametrically closer than in the WCP case. 

Comparing our theoretical results to typical measured 
current-voltage characteristics, we find good agreement 
with experimental results [17-21]. The excess-current 
characteristic reported in these experiments was pointed 



out early on by Campbell and Evetts[4], however, we 
are not aware of any 'microscopic' derivation of this ba- 
sic result. Unfortunately, a detailed comparison between 
theory and experiment is still not available today. Given 
a specific material, the defect structure is usually non- 
trivial and may include a variety of pin types. Further- 
more, the parameters characterizing the defects are dif- 
ficult to find. Experiments with superconductors where 
defects could be designed, tuned, and properly charac- 
terized would provide a great help and motivation in 
further developing the theory of pinning, particularly 
the crossover regime between strong and weak collective 
manifesting itself first in the small-velocity domain. 
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